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The relation between energy and density (known as the nuclear equation of state) plays a major 
role in a variety of nuclear and astrophysical systems. Spin and isospin asymmetries can have a 
dramatic impact on the equation of state and possibly alter its stability conditions. An example 
is the possible manifestation of ferromagnetic instabilities, which would indicate the existence, at 
a certain density, of a spin-polarized state with lower energy than the unpolarized one. This issue 
is being discussed extensively in the literature and the conclusions are presently very model depen- 
dent. We will report and discuss our recent progress in the study of spin-polarized neutron matter. 
The approach we take is microscopic and relativistic. The calculated neutron matter properties 
are derived from realistic nucleon-nucleon interactions. This makes it possible to understand the 
properties of the equation of state in terms of specific features of the nuclear force model. 
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1. INTRODUCTION 

The properties of dense and/or highly asymmetric nu- 
clear matter, where asymmetric may refer to isospin or 
spin asymmetries, are of great current interest in nuclear 
physics and astrophysics. This topic is broad-scoped 
since it reaches out to exotic systems on the nuclear chart 
as well as, on a dramatically different scale, exotic objects 
in the universe such as compact stars. 

In this paper, we investigate the bulk and single- 
particle properties of spin-polarized neutron matter. The 
study of the magnetic properties of dense matter is of 
considerable interest in conjunction with the physics of 
pulsars, which are believed to be rapidly rotating neu- 
tron stars with strong surface magnetic fields. The po- 
larizability of nuclear matter can have strong effects on 
neutrino diffusion and, in turn, variations of the neutrino 
mean free path due to changes in the magnetic suscepti- 
bility of neutron matter can impact the physics of super- 
novae and proton- neutron stars. 

The magnetic properties of neutron/nuclear matter 
have been studied extensively since a long time by 
many authors and with a variety of theoretical meth- 





theless, conclusions about the possibility of a phase tran- 
sition to a ferromagnetic state at some critical density 
are still contradictory. For instance, calculations based 
on Skyrme-type interactions [23j predict that such in- 
stabilities will occur with increasing density. In partic- 
ular, currently used Skyrme forces show a ferromagnetic 
transition for neutron matter at densities between l.lpo 
and 3.5/90 0|- On the other hand, more recent pre- 
dictions based on Monte Carlo simulations and the 
Brueckner-Hartree-Fock (BHF) approach with realistic 
nucleon-nucleon (NN) interactions exclude these 

instabilities, at least at densities up to several times nor- 
mal nuclear density. Similarly, no evidence of a transition 
to a ferromagnetic state was found in older calculations 
based on the Brueckner-Hartree-Fock approach with the 



Reid hard-core potential as well as a non-local separa- 
ble potentials Relativistic calculations based on ef- 
fective meson-nucleon Lagrangians [ISi] predict the ferro- 
magnetic transition to take place at several times nuclear 
matter density, with its onset being crucially determined 
by the inclusion of the isovector mesons. Clearly, the 
existence of such phase transition depends sensitively on 
the modeling of the spin-dependent part of the nuclear 
force and its behavior in the medium. Thus, this unset- 
tled issue goes to the very core of nuclear physics. 

Our calculation is microscopic and treats the nucleons 
relativistically. A parameter-free and internally consis- 
tent approach is important if we are to interpret our con- 
clusions in terms of the underlying nuclear force. This is 
precisely our focus, namely to understand the in- medium 
behavior of specific components of the nuclear force (in 
this case, the spin dependence). Different NN potentials 
can have comparable quality as seen from their global 
description of NN data and yet differ in specific features. 
Thus, it will be interesting to explore how, for a given 
many-body approach, predictions for spin-polarized neu- 
tron matter depend upon specific features of the NN po- 
tential. Second, it will be insightful to compare with pre- 
dictions based on a realistic NN potential and the BHF 
method [3], especially at the higher densities, where the 
repulsive Dirac effect can have a dramatic impact on the 
short-range nature of the force. 

This work is organized in the following way: after the 
introductory notes in this section, we briefly review our 
theoretical framework (section 2); our results are pre- 
sented and discussed in section 3; we conclude in section 
4 with a short summary and outlook. 



2. BRIEF DESCRIPTION OF THE 
CALCULATION 

The starting point of any microscopic calculation of nu- 
clear structure or reactions is a realistic free-space NN in- 
teraction. A realistic and quantitative model for the nu- 
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clear force with reasonable theoretical foundations is the 
one-boson-exchange (OBE) model [13] ■ Unless otherwise 
specified, our standard framework consists of the Bonn 
B potential together with the Dirac-Brueckner-Hartree- 
Fock (DBHF) approach to nuclear matter. A detailed 
description of our application of the DBHF method to 
nuclear, neutron, and asymmetric matter can be found 
in our earlier works ^32, 33] . 

Similarly to what we have done to describe isospin 
asymmetries of nuclear matter, the single-particle poten- 
tial is the solution of a set of coupled equations 

Uu = Uud + Uuu (1) 

Ud = Udu + Um (2) 



where u and d refer to "up" and "down" polarizations, 
respectively, and where each Uuui term contains the ap- 
propriate (spin-dependent) part of the interaction, Gaa' ■ 
More specifically. 



Ua{p)^ <a,a'\G{p,q)\o,a' >, (3) 

u'=u,d qKk'J 



where the second summation indicates integration over 
the two Fermi seas of spin- up and spin-down neutrons, 
and 



<a,a'\G{p,q)\n,a' > ^ ^ <K;]^cj'\S{a + a') ><]^a;]^a'\S{a + a') > 

L,L',S.,.J,M,Ml 

X < LML;S{a + a')\JM >< L' Ml; S{a + a')\JM > 

X i'''~''Yl,^j^j(Kel)YL.MAKel)<LSJ\G{kreuK,,^)\L'SJ> 



(4) 
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The notation < jimi; j2'7i2] jsm-a > is used for the 
Clebsh-Gordon coefficients. Clearly, the need to sepa- 
rate the interaction by individual spin components brings 
along angular dependence, with the result that the single- 
particle potential depends also on the direction of the 
momentum. Notice that the G-matrix equation is solved 
using partial wave decomposition and the matrix ele- 
ments are then summed as in Eq. (4) to provide the 
new matrix elements in the uncoupled-spin representa- 
tion needed for Eq. (3). The three-dimensional integra- 
tion in Eq. (3) is performed in terms of the spherical co- 
ordinates of g, (g, 9q, (t)q), with the final result depending 
upon both magnitude and direction of p. On the other 
hand, the scattering equation is solved using relative and 
center-of-mass coordinates, krei and Kc.m.- These are 
easily related to the momenta of the two particles in the 
nuclear matter rest frame through the standard defini- 
tions Kc.m. — P + Q, and krei = (The latter displays 
the dependence of the argument of the spherical harmon- 
ics upon p and q.) 

Solving the G-matrix equation requires knowledge 
of the single-particle potential, which in turn requires 
knowledge of the interaction. Hence, Eqs. (1-2) together 
with the G-matrix equation constitute a self-consistency 
problem, which is handled, technically, exactly the same 
way as previously done for the case of isospin asymmetry 
[3l| . The Pauli operator for scattering of two particles 
with unequal Fermi momenta, contained in the kernel of 
the G-matrix equation, is also defined in perfect analogy 
with the isospin-asymmetric one [sij ]. 



Qaa' (Pi 9j kp, kp ) 



1 \ip> kp and q > kp 
otherwise. 



(5) 



Notice that, although a full three-dimensional integration 
is performed in Eq. (3), the usual angle-average proce- 



dure is applied to the Pauli operator (when expressed in 
terms of krei and Kc.m.) and to the two-particle propa- 
gator in the kernel of the G-matrix equation. 

Once a self-consistent solution is obtained for the 
single-particle spectrum, the average potential energy for 
particles with spin polarization a is obtained as 
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2 (27r)3 p. 



p=0 



Ua{p)p^dpd^p 



J 4,^=0 



The average potential energy per particle is then 

y_ P« <Uu> +Pd < Ud > 



(6) 



(7) 



The kinetic energy (or, rather, the free-particle opera- 
tor in the Dirac equation when using the DBHF frame- 
work), is also averaged over magnitude and direction of 
the momentum. In particular, we calculate the average 
free-particle energy for spin-up(down) neutrons as 



<T„ > = 



fdnf{mue)) 
jdn 



(8) 



where T is the average over the magnitude of the mo- 
mentum. Notice that the angular dependence comes 
in through the effective masses, which, being part of 
the parametrization of the single-particle potential, are 
themselves direction dependent (and of course different 
for spin- up or spin-down neutrons). 
Finally 

< T >^ Pu<Tu> +Pd <Td> 
P 

and the average energy per neutron is 

e=<U> + <T> (10) 
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FIG. 1: (Color online) Angular dependence of the single- 
particle potential for spin-up and spin-down neutrons at fixed 
spin asymmetry and Fermi momentum and for different val- 
ues of the neutron momentum. The momenta are in units of 
fm~^. The angle is defined relative to the polarization axis. 



As in the case of isospin asymmetry, it can be expected 
that the dependence of the average energy per particle 
upon the degree of polarization [18] will follow the law 



e(p,/3) = e(p,/3-0) + 5(p)/32 



(11) 



where (3 is the spin asymmetry, defined hy P — — — —. 
A negative value of S{p) would signify that a polarized 
system is more stable than unpolarized neutron matter. 
From the energy shift. 



^(p) = e(p,/3=l)-e(p,/3 = 0). 



(12) 



the magnetic susceptibility can be easily calculated. If 
the parabolic dependence is assumed, then one can write 
the magnetic susceptibility as 



X 



25(p)' 



(13) 



where fj, is the neutron magnetic moment. The magnetic 
susceptibility in often expressed in units of xf, the mag- 
netic susceptibility of a free Fermi gas. 



XF 



(14) 



where k p denotes the average Fermi momentum which is 
related to the total density by 

kp = (3^2^)1/^ (15) 

The Fermi momenta for up and down neutrons are 



k'^p = kp{i + i3y/^ 
kf ^ kp{i ~ i3y/\ 



(16) 



FIG. 2: (Color online) As in the previous figure but for a 
larger value of the spin asymmetry. 



For the most general case, it will be necessary to com- 
bine isospin and spin asymmetry. With twice as many 
degrees of freedom, the coupled self-consistency problem 
schematically displayed in Eqs.(l-2) is numerically more 
involved but straightforward. This is left to a future 
work. 



3. RESULTS AND DISCUSSION 

We begin by showing the angular and momentum de- 
pendence of the single-neutron potential, see Figs. 1-2. 
The angular dependence is rather mild, especially at the 
lowest momenta. As can be reasonably expected, it be- 
comes stronger at larger values of the asymmetry, see 
Fig. 2. In Fig. 3, the asymmetry dependence is displayed 
for fixed density and momentum (here the angular de- 
pendence is averaged out). As the density of u particles 
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FIG. 3: (Color online) Asymmetry dependence of the single- 
particle potential for spin-up and spin-down neutrons at fixed 
density and momentum. The angular dependence is inte- 
grated out. 
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FIG. 4: (Color online) Asymmetry dependence of the effective 
masses for upward and downward polarized neutrons under 
the same conditions as in Fig. 3. 
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FIG. 5: (Color online) Average '^potential, kinetic, and to- 
tal energy per particle at various densities as a function of 
the spin asymmetry. Predictions obtained with our standard 
DBHF calculation. 
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goes up, the total density remaining constant, the most 
hkely kind of interaction for u neutrons is of the uu type. 
Similarly, the largest contribution to the d-particle po- 
tential is of the du type, see Eqs. (1-2), with the latter 
being apparently more attractive, as can be inferred by 
the spin splitting of the potential shown in Fig. 3. Be- 
fore we move on to discuss this issue in greater details, 
we also show the effective masses of u and d neutrons, 
see Fig. 4, and observe that they display a qualitatively 
similar behavior as the one of the corresponding single- 
particle potentials. 

The average energy per particle at various densities 
and as a function of the asymmetry parameter is shown 
in the third frame of Fig. 5. The first two frames display 
the contribution from the average potential energy and 
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FIG. 6: (Color online) Same legerid as for the previous figure, 
but the predictions are obtained with the BHF calculation. 
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FIG. 7: (Color online) Energy difference between the polar- 
ized and the unpolarized states corresponding to Fig. 5 and 
6. 



the average kinetic energy, respectively. The parabolic 
dependence on f3, or linear on is obviously verified. 
In Fig. 6 we show the corresponding predictions obtained 
with the conventional Brueckner-Harteee-Fock approach. 
This comparison may be quite insightful, as we further 
discuss next. We notice that the Dirac energies are over- 
all more repulsive, but the parabolas predicted with the 
BHF prescription appear to become steeper, relative to 
each other, as density grows. The energy difference be- 
tween the totally polarized state and the unpolarized one 
for both the relativistic and the non-relativistic calcula- 
tion is shown in Fig. 7. Although initially higher, the 
growth of the DBHF curve shows a tendency to slow 
down and the two sets of predictions cross over just above 

3/30 . 

Before leaving this detour into the non-relativistic 
model, we observe that the predictions shown in Fig. 6 
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FIG. 8: (Color online) Neutron effective masses used in the 
DBHF calculations of the EOS. The angular dependence is 
averaged out. 
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FIG. 9: (Color onhne) Average energy per particle at densi- 
ties equal to 0.5,1,2,3,5,7,9, and 10 times po (from lowest to 
highest curve). Predictions obtained with our DBHF calcu- 
lation. 



concerned, we would expect them to be more pronounced 
for nuclear matter than for pure neutron matter, since 
the largest variations among modern realistic potentials 
are typically found in the strength of the tensor force, 
which is stronger in T=0 partial waves (obviously absent 
in the nn system). This point will be explored in a later 
investigation. 

In the remainder of this paper, we will focus on the 
DBHF model, which is our standard operational ap- 
proach. To further explore the possibility of a ferro- 
magnetic transition, we have extended the DBHF cal- 
culation to densities as high as lOpo- The same method 
as described in Ref. [s^l is applied to obtain the energy 
per particle where a self-consistent solution cannot be 
obtained (see Section HI of Ref. jsl] for details). The 
(angle-averaged) neutron effective masses for both the 
unpolarized and the fully polarized case are shown in 
Fig. 8 as a function of density. 

DBHF predictions for the average energy per particle 
are shown in Fig. 9 at densities ranging from p=0.5po to 
lOpo- What we observe is best seen through the spin- 
symmetry energy, which we calculate from Eq. (12) and 
show in Fig. 10. We see that at high density the energy 
shift between polarized and unpolarized matter continues 
to grow, but at a smaller rate, and eventually appear to 
saturate. Similar observations already made in conjunc- 
tion with isospin asymmetry were explained in terms of 
stronger short-range repulsion in the Dirac model [ssj . 
It must be kept in mind that some large contributions, 
such as the one from the ^Sq state, are not allowed in 
the fully polarized case. Now, if such contributions (typ- 
ically attractive at normal densities) become more and 
more repulsive with density (due to the increasing impor- 
tance of short-range repulsive effects), their absence will 
amount to less repulsive energies at high density. On the 
other hand, if large and attractive singlet partial waves 
remain attractive up to high densities, their suppression 
(demanded in the totally polarized case) will effectively 
amount to increased repulsion. 

In conclusion, although the curvature of the spin- 
symmetry energy may suggest that ferromagnetic insta- 
bilities are in principle possible within the Dirac model, 
inspection of Fig. 10 reveals that such transition does 



are reasonably consistent with those from previous stud- 
ies which used the Brueckner-Hartree-Fock approach and 
the Nijmegen II and Reid93 NN potentials In fact, 
comparison with that work allows us to make some useful 
observations concerning the choice of a particular NN po- 
tential, for a similar many-body approach (in this case, 
BHF). We must keep in mind that off-shell differences 
exist among NN potentials (even if nearly equivalent in 
their fit of NN scattering data) and those will impact the 
G-matrix (which, unlike the T-matrix, is not constrained 
by the two-body data). Furthermore, off-shell differences 
will have a larger impact at high Fermi momenta, where 
the higher momentum components of the NN potential, 
(usually also the most model dependent), play a larger 
role in the calculation. Accordingly, the best agreement 
between our BHF predictions and those of Ref. 18] is 
seen at low to moderate densities. Furthermore, as far as 
differences based on the choice of the NN potential are 
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FIG. 10: (Color online) Density dependence of the spin sym- 
metry energy obtained with the DBHF model. 
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CONCLUSIONS 
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FIG. 11: (Color online) Density dependence of the ratio 
Xf/x- Predictions are obtained with the DBHF model. 



not take place at least up to lOpo- Clearly it would not 
be appropriate to explore even higher densities without 
additional considerations, such as transition to a quark 
phase. In fact, even on the high side of the densities 
considered here, softening of the equation of state from 
additional degrees of freedom not included in the present 
model may be necessary in order to draw a more definite 
conclusion. 

Finally, in Fig. 11 we show the ratio xf/Xj whose be- 
havior is directly related to the spin-symmetry energy, 
see Eq. (13). Clearly, similar observations apply to both 
Fig. 11 and Fig. 10. (The magnetic susceptibility would 
show an infinite discontinuity, corresponding to a sign 



We have calculated bulk and single-particle proper- 
ties of spin-polarized neutron matter. The EOSs we ob- 
tain with the DBHF model are generally rather repulsive 
at the larger densities. The energy of the unpolarized 
system (where all nn partial waves are allowed), grows 
rapidly at high density with the result that the energy 
difference between totally polarized and unpolarized neu- 
tron matter tends to slow down with density. This may 
be interpreted as a precursor of spin-separation instabil- 
ities, although no such transition is actually seen up to 
lOpo- Our analysis allowed us to locate the origin of this 
behavior in the contributions to the energy from specific 
partial waves and their behavior in the medium, partic- 
ularly increased repulsion in the singlet states. 

In future work, the impact of further extensions will 
be considered, such as: examining the effects of con- 
tributions that soften the EOS (especially at high den- 
sity) ; extending our framework to incorporate both spin- 
and isospin-asymmetries; examining the temperature 
dependence of our observations for spin- and isospin- 
asymmetries of neutron and nuclear matter. 
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